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Abstract 

We construct two new classes of models for double beta decay, each of which leads to an 
electron spectrum which differs from the decays which are usually considered. One of the 
classes has a spectrum which has not been considered to date, and which is softer than 
the usual two-neutrino decay of the Standard Model. We construct illustrative models to 
show how other phenomenological bounds can be accommodated. We typically find that, 
although all other bounds can be satisfied, the predicted double-beta decay rate only in one 
class of these models is at best large enough to be just detectable in current experiments. 
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1. Introduction 



The purpose of this paper is to present two classes of double-beta decay, which 
have not been hitherto considered by workers in the field, and are potentially experi- 
mentally distinguishable from the presently-known kinds of decays. Besides presenting a 
general discussion of the phenomenology oi (3(3 decay in the kinds of theories which we 
consider, we also explore the other kinds of experiments to which these models might be 
expected to contribute. In so doing we fill in the last missing items in the classification of 
(3(3 decay that was presented in ref. [1]. 

First, a brief motivation for this calculation. 

The past ten years have seen a great deal of effort go into the experimental search 
for (3(3 decay, an extremely rare process in which two neutrons simultaneously decay into 
protons with the accompanied emission of two electrons. This process — when it oc- 
curs together with the emission of two antineutrinos {(3(32i,) — is predicted to occur at 
second-order in the charged-current weak interactions in the Standard Model (SM). The 
experimental effort has borne fruit in recent years, with the unambiguous observation of 
this decay in a number of different nuclear isotopes [2]. 

Of course, part of the interest in these experiments was originally motivated by the 
possibility of observing other decays, besides the SM process. There are several other 
types of decay modes to which these experiments are sensitive, such as the lepton-number 
violating neutrinoless process (/?/?oi/), or the decay {(3(3^p) into two electrons plus a neutral 
scalar, (p. Both of these reactions have been predicted by a number of plausible alternatives 
to the SM, starting with the simple Gelmini-Roncadelli (GR) model of lepton-number 
breaking [3], [4]. 

The experimental quantity that is used to distinguish such exotic decays from the run- 
of-the-mill SM events is the shape of the decay rate as a function of the energies, Si and 
£2, of the two emitted electrons. This shape is simplest for (3(3qi, decay, since in this case 
the rate is zero unless the sum of the two electron energies, 6 = 61+62, takes on a specific 
value, Q. Q represents the total energy that can be carried off by the two electrons, and is 
given in terms of the relevant masses — denoted respectively by M and M' for the parent 
and daughter nuclei — by Q = M — M'. In all of the decays of interest (5 — 2 MeV in 
size. Other decays (like /?/?2z/ and /3/?<p), which involve other relativistic decay products in 
addition to the final electrons, instead predict a continuous decay distribution throughout 
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the entire interval 2me < s < Q. 

The spectrum that is predicted by the various types of decays turns out, quite gener- 
ally, to take a particularly simple form that is completely characterized by a single integer, 
or 'spectral index', n. Explicitly: 

_iL_ = c{Q - ei - e^r \plelF{e^)] \p2e2F{e2)] , (1) 

where C is independent of Si and 62, and Pi = \pi\, for i = 1,2, represents the magnitude 
of the three- momentum of the corresponding electron. The quantity F{si) is the Fermi 
function which describes the spectral distortion due to the electric charge of the nucleus. 
The spectral shape which follows from eq. (1) is shown, for various choices of n, in Fig. (1). 

It is ultimately the small size of the energy, Q ~ 2 MeV, in comparison with the typical 
momenta, pp ~ 60 MeV, of the nucleons in the decaying nucleus that ensures that the 
electron spectrum takes such a simple form — i.e. one that is characterized simply by the 
integer n. This is because the small ratio of these two scales justifies keeping only the lowest 
term in an expansion of the decay rate in powers of the momenta of the final electrons 
(and/or scalars and neutrinos). The index, n, that is appropriate for any particular kind of 
decay is determined by the leading term in this expansion of the relevant decay rate. For 
example, for the SM process, PP2U, the phase space volume element for the final neutrinos 
implies an index n{SM) = 5. For the scalar decay, of the GR model, on the other 
hand, the phase space of the lone scalar implies n{GR) = 1. The diff'erence between these 
predictions — c.f. Figure (1) — forms the basis for the experimental discrimination of these 
two models. 

Two developments have recently provoked a theoretical re-evaluation [5], [6], [7] of 
the kinds of new physics to which the current experiments can be sensitive. The first of 
these was the advent of precision measurements of the properties of the Z resonance at 
LEP [8], which has ruled out the original GR-type models for new physics. The other 
development has come from the experiments themselves, starting with the appearance 
of several indications of an excess of electrons just below the endpoint in some of the 
experimental data for several types of decays [9] . Although the evidence for the excess has 
since diminished [10], the conclusion reached by the theoretical re-evaluation still stands: 
viable theories which can produce observable decay are possible, and can be consistent 
with all other, non-/3/?, experimental limits, but only if the new physics has rather different 
properties than have previously been assumed based on experience with GR-type models. 
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One of the main surprises which arose from this theoretical re-examination has been 
the observation that new physics can be very generally divided into a small number of 
classes — depending on how it addresses a few basic issues — and that the experimental 
features that any model predicts for (3(3 decay is a robust indicator of the class to which it 
belongs [1]. (These issues and categories are summarized in Table I.-^ ) 
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Table I 

A list of alternatives for modelling double beta decay. 

As is argued in ref. [1], the observed absence of /?/?oi/ decays at the endpoint of the 
electron spectrum strongly constrains the classes of models for which lepton number is 
broken {i.e. classes I A - IE). In particular, it implies that all of the scalar-emitting decays 
in these classes are essentially indistinguishable from their counterparts in the lepton- 
number preserving classes. That is to say, so long as PPou decays are not observed, models 

^ The table presented here differs sUghtly from the table in ref. [1] . A trivial change is our introduction 
of new categories ID and IE, even though these are in practice indistinguishable from IIC. Also, the index 
listed here for category HE reflects the results of the present paper, and differs from the preliminary guess 
for this value which was given in ref. [1]. 
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in class IB and IC are at present experimentally indistinguishable from those in class IIB. 
Similarly, classes ID and IE cannot be distinguished from class IIC. 

Interestingly, all of the models which had been considered previously, as well as most 
of the more recent proposals, fall into only a few of the several possible categories — cases 
IB and IC of Table I. The other, unorthodox, categories remain essentially unexplored. 
The two exceptions to this statement are case IID, which was discovered and studied in 
ref. [5], and case IE, which contains a supersymmetric model [11], in which two scalars are 
emitted m. (3(3 decay: PP^ip. Unfortunately, this last model has since been ruled out by 
the measurements at LEP. Both of these new classes of models predict a /?/?<p spectrum 
having n = 3. The resulting shape is intermediate between the usual SM and GR spectra, 
and can be distinguished from these in current experiments [12], [13]. 

The purpose of the present paper is to explore in detail those classes of theories of 
Table I which have not been considered to date. Keeping in mind the indistinguishability 
of classes IB, IC and IIB, as well as ID, IE and IIC, we see there are two types of models 
which we must here consider: those of classes IIC and HE, both of which involve decay 
accompanied by the emission of two, rather than just one, scalar particles. Two scalars 
must be emitted because of lepton-number conservation, and the charge assignments of 
the light particles. Models in these two categories differ from one another according to 
whether or not the light scalar particle that is emitted is a Goldstone boson. 

We first examine models of category IIC to show that these can be altered to evade 
detection at LEP, and we then work through the completely new category, HE, which 
predicts a spectral index, n = 7, that is completely different from all previous cases. Since 
this index is larger than that for the SM decay — n{SM) = 5 — the resulting spectral 
shape is softer than the presently-observed f3f32v spectra. This kind of softer spectrum is 
much more diflBcult to distinguish from the experimental background, and so we expect 
the experimental limit for the branching ratio into the new n = 7 decay mode to be 
comparatively weak. 

On the theoretical side, we find that the nuclear form factors which appear in all of 
the new /?/?<p<p decays are well understood, since they are the same as those that appear 
in the more well-known decay modes, such as for /?/3oi/ and /3/?<p in the GR model. We 
find that using these matrix elements to estimate the potential size of the /3/?<p<p decay rate 
typically gives a result which can be close to the present experimental limit for the n = 3 
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decays, provided we relax the conditions that are cast on the relevant masses by big bang 
nucleosynthesis constraints. The estimated rate for the n = 7 decays appears to be too 
small to be detectable. 

This paper has the following organization. The next section computes the /3/3<^(^ decay 
rate for the two classes of decays (IIC and HE) . We do so to verify their spectral index — 
which is their experimental signature — and to determine how the resulting rates depend 
on the couplings in underlying models. For class IIC (having n = 3) we consider two types 
of decay mechanisms: those which proceed due to the exchange of new sterile neutrinos, 
and those which proceed through the interactions of new scalar particles. Representative 
models for decays in class IIC are then constructed in sections 3, and those for class HE 
in section 4. Section 5 then summarizes the most important phenomenological bounds for 
these models. Our conclusions are finally summarized in section 6. 

2. Two-Scalar pp Decays 

Before constructing models in detail, we first pause here to record expressions for the 
decay rate into two electrons plus two scalars, T{PP^^). We do so in order to display how 
the PPiptp rate depends on the various masses and couplings of the theory. We may then use 
this dependence to see what kinds of new particles can be used to generate an observable 
signal in experiments. We defer the discussion of any further phenomenology, which 
requires the construction of an explicit model, to the following sections. 

We consider here three kinds of scenarios. First we consider two types of models for 
which (3(3^ip decay is mediated by the admixture of ordinary neutrinos with new sterile 
fermions. The two kinds of such decays which we consider are those for which the scalar 
emission is either derivatively suppressed {i.e. case HE, with n = 7), or not (case IIC with 
n = 3). Finally we consider the alternative for which /3/?<py, decay occurs because of the 
interactions of new types of scalar particles. Such scalar-mediated decays can also arise 
with or without derivative suppression for the final scalar emission (classes HE or IIC), 
but since the decay rate for the class HE models are generally very small, we do not pursue 
them in detail for the scalar-mediated case. 

2.1) Fermion-Mediated n = 3 Decays 

Consider first a theory of neutrinos, Vi and Na, which respectively carry lepton number 
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Le{^i) = +1 cind Lf,{Na) = 0, and which are coupled to a scalar, 0, which has lepton num- 
ber l/e(0) = +1. The most general renormalizable and Lg-conserving Yukawa couplings 
involving these fields is: 



^ynk ^ -^i {AiajL + Bia-fn) Na ((> + h.C, (2) 

where Aia and Bia represent arbitrary Yukawa-coupling matrices, and 7^ and 7h denote 
the usual projections onto left- and right-handed spinors. We use majorana spinors here 
to represent the neutrinos, so the conjugate which appears in eq. (2) is Vi = uj C~^, where 
C is the charge-conjugation matrix. Suppose also that a number of right-handed neutrinos 
are included in the theory, in order to form generic lepton number conserving masses, rrtj/. , 
for the Le = 1 states, and that the Le = neutrinos, Na, have generic majorana masses. 

As long as any of the participate in the charged-current weak interactions, the 
couplings of eq. (2) generically produce f3(3^p^p decay due to the Feynman graph of Fig. (2). 
(Two scalars must be emitted in this decay due to conservation of Lg.) Even though Fig. (2) 
arises at tree level, the four-momentum of the exchanged neutrino is not determined by 
energy and momentum conservation. This is because the amplitude for the decay of the 
two initial neutrons must be convoluted with the wavefunction of the initial nucleus. The 
important momentum scale in the integration over the exchanged neutrino is therefore set 
by the scale of momenta to which the nucleons within the initial and final nuclei have access. 
Since this scale is much larger than the energy, Q, that is available for the final electrons 
and scalars, it is a good approximation to neglect the outgoing scalar and electron energies 
in the decay amplitude. Evaluating the graph in this approximation gives the following 
expression for the decay rate: 

dV{(3(3^^) ^^^^^^ \A{(5(3^^f {Q-ei-e2fY[pkekF{ek) dsk, (3) 

where Gp is Fermi's constant, and 9c is the Cabbibo angle. The quantity A{(3I3^^) repre- 
sents the integral: 



I] a 



APP^^) [^^^2) J ^27r)^{e^ + ml- ie){P + ml.- ie){P + ml -iey 
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where Vei/j is the Kobayashi-Maskawa-type matrix for the leptonic charged current. The 
factor, Nijai in the numerator of eq. (4) denotes the foUowing expression: 




(5) 



Finally, the tensor W^i^ is an independent function of the space- and time-like parts, 
£^ and |1|, of £^ in the rest frame of the decaying nucleus, which encodes the nuclear matrix 
element that is relevant to the decay [5] : 



Here = 2ll7^7j,(i denotes the appropriate hadronic charged current. The prefactor in- 
volving the energy-to-mass ratio — E/M and E' /M' , for the parent and daughter nuclei 
respectively — is required in order to ensure that Wf^iy transforms under Lorentz transfor- 
mations as a tensor. The factor of (27r)^ is purely conventional. 

There are several features in these expressions that are noteworthy. 

1. Comparison of eqs. (3) and (1) shows that the spectral index for this decay is n = 3. 
The additional two powers of {Q — e), in comparison with the GR model, arise here 
simply due to the additional phase space volume of the second scalar. 

2. Eq. (4) depends on precisely the same combination of nuclear form factors, W^^/j, = 
Wp — Wgti as do f3f32v and f3f3ov decays, as well as f3f3^p decay in GR-type models 
[5] . (For future use we express ^j, in terms of the Fermi and Gamow-Teller matrix 
elements, Wp = Wqq and Wgt = Yli=i^ii- R-^f- [5] gives Wp- and explicit 
expressions in terms of the nuclear matrix elements of various neutrino potentials.) 
Since these particular combinations of nuclear matrix elements have been well studied 
within the context of these other processes [14], [15], [16], there is relatively little 
uncertainty in the estimate for the total rate for this type of (3(3ipip decay. 

3. The differential decay distribution for this decay as a function of the opening angle, 9, 
between the two electrons is exactly the same as it is for (3(32u, P(3oi/ and /3/3<^ decays 
(of both the GR type and the new (3(3^ decays of type IID). It is given explicitly 
(neglecting the mutual repulsion of the outgoing electrons) by 



d^'x (iV'|T* [J^{x)MO)] \N) e 



(6) 



1 dT _1 
f dcos^~ 2 



(1 — V1V2 cos 9) , 



(7) 
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where Vi = Pi/si, for i = 1,2, is the speed of the corresponding electron in the nuclear 
rest frame. 



2.2) Fermion-Mediated n — 7 Decays 

We next consider the case of double scalar emission, but in the special case where 
the amplitude for emitting the scalars is suppressed by the scalar energy. This is the case 
when the scalars have purely derivative couplings, such as when the scalar is a Goldstone 
boson which carries an unbroken lepton number [5] . For such models the expression given 
in eq. (4) for the matrix element. A, vanishes identically, and so it is necessary to work to 
higher order in the external scalar momenta in order to obtain a nonvanishing rate. Since 
it turns out that the first nonvanishing contribution to the amplitude arises at quadratic 
order in the scalar momenta, the calculation is quite cumbersome when performed directly 
with the couplings of eq. (2). 

A better way to proceed is to perform a change of variables so that the derivative 
coupling nature of the Goldstone bosons is manifest from the outset. Once this has been 
done, the trilinear coupling to neutrinos of a Goldstone boson carrying Lg = 1, becomes: 



where the coefficients Xia and Yia are computable [5] in terms of Aia and Bia of eq. (2), as 
well as the components of the broken-symmetry generator, Q, for which is the Goldstone 
boson. 

The contribution to PPip,p decay in these models arises from the same Feynman graph 
as before, Fig. (2), but with the interaction of eq. (8) at each vertex. Evaluation of this 
graph gives a fairly complicated expression for the decay rate, but this expression simplifies 
greatly in the special case of purely left-handed couplings, for which Yia = 0. Since this 
case is sufficient to analyse the models we consider in later sections, we just quote here the 
decay rate in the limit Via — 0: 




(8) 



{G^ cos 6>c)^ 
8(27r)5 



MPP^^) (Q-ei- 62)' n Pk^kFisk) dsk, (9) 
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where the new quantity, A{PP^p^), represents the integral: 

MPS )=(^)'Tf^ v;..Ve.,M> 

■^{PPvv) \^io5^2^ Z^J (27T)'(e^ + ml-i€)(P+mi.-i€)(P+ml^-i€y ' 

Apart from its normalization, eq. (10), differs from eq. (4), obtained previously, only by 
the replacement Mija J^ija, where: 

Afija={-£^)XiaXjam^^. (11) 

Once again there are several properties of the above expressions that bear emphasis. 

1. The spectral index for this decay is n = 7. The four additional powers of {Q — e), 
in comparison with eq. (3), arise because each scalar-emission vertex of eq. (8) is 
explicitly proportional to the scalar four-momentum. 

2. The decay rate depends only on the form factor combination W^^ = Wp- — Wgt, just 
as for the non-derivatively coupled calculation. (This is by contrast with the situation 
for 13(3^ decay, where the matrix elements for the derivatively-suppressed decay — 
class IID — differ from those for the emission of a generic scalar — class IIB.) 

3. The differential decay distribution as a function of the electron opening angle remains 
exactly the same as for all of the other decays, as in eq. (7). 

2.3) Scalar-Mediated n = 3 Decays 

The previous calculations produced similar expressions for the decay amplitudes — 
c.f. eqs. (4) and (10) — largely because the same Feynman graph is responsible for it in 
both classes of models. In this section we explore an alternative class of models for which 

(3(3ipip decay proceeds via a different process: the exchange of an exotic set of scalars. 

Consider therefore supplementing the SM by a collection of new scalars, which we 
classify using their lepton and electric-charge quantum numbers, Lg and Q. Besides the 
very light (Q, Lg) = (0, -|-1) particle, (f, which is emitted in the decay, we imagine also 
adding two other new classes of fields, Xi ^^d A^, which respectively have the {Q,Le) 
assignments: (-|-1,— 1) and {—2, +2). These quantum numbers permit the following types 
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of Lg-invariant trilinear-scalar, and Yukawa couplings: 




-e(/i27^ + /i«7«)e'^Aa + h.c., 



(12) 



(13) 



in addition to the charged-current coupling 



(14) 



Here e'^ denotes the Dirac conjugate of the electron field, and A' parameterizes the strength 
of the scalar charged-current interaction relative to the usual SUj^{2) gauge coupling, g. 
For the sake of keeping our final expressions simple, we assume from here on that the 
index 'a' takes only one value, so that there is only one field. A, having the quantum 
numbers {Q,Le) = (—2,-1-2). We therefore omit this index from the coupling parameters 
as these are defined above: i.e. — — > /i*-', and h'^ j^ — > h^^ji. We also denote the mass 
eigenvalues for Xi and A respectively by Mj and Ma- 

With these couplings and interactions decay is mediated by the Feynman graph 

of Fig. (3). Evaluation of this graph gives the following decay rate: 



dr.ipp^^) = ^ I \As{f3f3^^)\' {Q-si- e^fUpkSkFiek) dsk, (15) 



2 



k=l 



where the scalar- mediated decay amplitude, As{P(3ipip), is: 




(27r)4 (£2 + M2 - ie) {P + Mf - 



ie) 



(16) 



In this expression h = \^\h^J^~+~\h^ ■, and 




(17) 



with ^ = 2Re [hLh%] and me denoting the electron mass. 
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These expressions differ from those for the fermion-mediated decays in a number of 
interesting ways. 

1. The spectral index for this decay is determined by phase space, and is n = 3. This is 
because the assumed coupHngs are not derivatively suppressed. A class of similar 
n = 7 decays can be obtained by imposing the additional requirement of derivative 
coupling. We do not consider such decays explicitly here. 

2. Once more the decay rate depends only on the familiar form factors, Wp and Wct, 
although this time they do not arise in the particular linear combination W^^ = 

Wp — Wqt- 

3. If the A-electron coupling has a right-handed component, 7^ 0, then the differential 
decay distribution as a function of the electron opening angle differs from that of all 
of the other decays we have considered. Assuming only one species of A, its explicit 
expression is. 



so the deviation of the decay distribution from that of eq. (7) is controlled by the 
difference between R and unity. 

We next turn to the construction of some underlying models for each of these types 



We first consider models for which the dominant contribution to (3(3 decay is through 
the emission of two scalars, but for which the emission of these scalars is not suppressed 
by powers of the scalar momentum, i.e. models of type IIC. As was mentioned earlier, 
and as is discussed in some detail in refs. [1] and [5], this also includes those models for 
which the light scalar is the Goldstone boson for lepton number (class IE), even though 
the couplings of such a scalar can always be put into the derivative form of eq. (8).^ 

In the following sections we consider a few representative models which produce this 
kind of decay. 

Briefly, with derivatively-coupled variables the dominant contribution to /3/3 decay for small scalar 
momentum comes from graphs in which the scalar is emitted from the electron line, for which an infrared 
singularity compensates for the derivative coupling. 




(18) 



of decays. 



3. Model-Building: The Case of Spectral Index n = 3. 
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3.1) Decays Mediated by Sterile Neutrinos 

A mechanism for producing PP^^ decay, without running into conflict with the LEP 
bounds, is to require the hght scalar to be an electroweak singlet, which couples to the usual 
electroweak neutrinos through their mixing with various species of sterile neutrinos. This 
results in a variation of the old singlet- majoron model [17]. In this way an unacceptable 
contribution to the invisible Z width can be avoided. A simple realization of this idea 
requires the introduction of three species of electroweak-singlet, left-handed neutrinos: N± 
and A'^o, together with a single scalar field, (p, carrying lepton number +1. The subscript 
i — 0,± on Ni indicates the lepton number of the corresponding left-handed state. 

Ignoring, for simplicity, all mixings with and i/r, we are led to the following, most 
general, lepton-number-invariant mass and Yukawa terms involving the new particles: 

Cm = -M(Ar+7,Ar_) --(AToT.Ao) +h.c., ^^^^ 
Ly = -A(]V_7,L) H-g_ (N_^,No) ^ - g+ (N+l,No) 0* + h.c. 

Here L = {^^^ and H = (^o) represent the usual SM lepton and Higgs doublets. (Where 
necessary the conjugate Higgs doublet is denoted by ^ = ia2H* .) Finally, recall that our 
use of majorana spinors implies Ni = N[C~^. 

Although, in general, both and (p could acquire vacuum expectation values (vef's): 
if'^) = V and (0) = in practice it is a good approximation to work in the limit u ^ 0. 
This is because the absence of an observable (3(3oi, decay signal implies the combination 
giu must be small. On the other hand, the couplings gi themselves cannot be too small or 
else the PP^p^p decay of interest here will itself be unobservable. As a result the vev u must 
itself be constrained to be negligible in comparison to the other masses of the problem. We 
therefore neglect u in all of what follows.^ We nevertheless assume that the scalar, 0, is 
light enough to permit PPypyp decay to occur. This can be done by tuning the mass in the 
potential to be small, or by permitting u ^ but with u tuned to be small enough to not 
produce too much PPoi^ decay. Either choice involves a certain amount of arbitrariness. 

The neutrino mass spectrum for this model greatly simplifies in the u = limit. There 
is one Lg — mass eigenstate, Nq, with mass m; there is a massless state, u'^ — UeCo — N^sg; 

o 

Notice that, unlike the choice ^=0, a very small but nonzero value for u is quite difficult to make 
natural, even in the technical sense [5], [18]. 
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and there is a massive, Dirac neutrino, Vg = IlN'^ + 7j^A^_, with A^^ = N^cq + i/gSg and 
having mass Mg = y/M^ + X^v^. Here ce and se denote the cosine and sine of the mixing 
angle, 6, which is defined by: tan^ = Xv/M. 

The Yukawa couphngs of these mass eigenstates take the form of eq. (2), with 

^^No = 0, A„^No =9-, -Si/^jvo = -g+se, B„^tfo = g+ce- (20) 
Using these expressions in the general result, eq. (4), then gives: 

^IPP.^.,; - |^3^2 ) '""'s J ^2^^4 (£2 _ (£2 + M| - ie)2 + rn?- ie) ' ^ ^ 

This expression implies that the I3(3^p^p decay rate is maximized for large couplings 
and mixing angles, and for sterile-neutrino masses in the vicinity of the nuclear-physics 
scale, ~ (10 — 100) MeV, which defines the important integration region in eq. (21). In 
order to see how big this rate can get we therefore have to see how close we can get to 
this optimal range using phenomenologically acceptable values for the various parameters. 
The relevant constraints are summarized in section 5, and they suggest that the amplitude 
can be made biggest for the following values of masses and mixing angles: 

Ms ~ 350 MeV; m ~ 1 MeV; ^ ~ 7.5 x 10"^; g+ ~ 0.2 (22) 

or, provided we relax the bound on m coming from nucleosynthesis (as explained in section 
5) 

Ms ~ m ~ 350 MeV; ^ ~ 7.5 x 10"^; g+ ~ 1 (23) 

Interestingly, the coupling constant g- need not be large. Notice also that the smallness 
of 9 implies Aw < Mg fa M. 

Using these values, and recalling that the integration momentum satisfies ^ < 60 MeV <S 
Ms, we find eq. (21), to be given approximately by 

/ 2 \^ f d^£ 
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For historical reasons, experimentalists conventionally state their results concerning 
the nonobservation of a scalar-emitting decay in terms of the old Gelmini-Roncadelli type 
j3j3^f, decay. The quantitative limit is expressed in terms of an upper bound to a dimension- 
less effective coupling constant, ^eff, which parametrizes the electron- neutrino/majoron 
interaction in the GR model: 

>Cphen = ^^i^e75J^e (25) 

This interaction gives a fifi^ decay rate which can be obtained from eq. (3) by performing 
the replacement: A{I3I3^^) {Q — Si — £2)^ A{GR) {Q — Si — £2), with 

Comparing this prediction with the results of current experiments gives an upper bound 

^eff < 10-^ [13]. 

To obtain a rough estimate of the sensitivity of current experiments to the PP^ip signal 
predicted by the model, we compare the PP^p^ decay rate with the /?/?<p rate which would 
be just detectable for the GR model, given the current experimental limit ^gfi ^ 10~^. For 
this kind of comparison the precision of a detailed matrix element analysis is unnecessary, 
so we follow Refs. [18] and [19] and make the following two approximations, which permit 
an analytical calculation of the various rates. 

1 To perform the integral over £, we approximate the nuclear matrix elements by 
replacing the form factors, Wp and Wot by step functions in energy and momen- 
tum: Wi ^ w^Q{Ep. — Io)Q{pp — |1|). decay rates are reproduced if we take 
\w^ — w'^r^l ~ 4 MeV""*^, and the present experimental sensitivity to an electron- 
neutrino majorana mass in PPoi, decay requires a nuclear Fermi momentum of ~ 60 
MeV. The Fermi energy then becomes Ep = p^/2mjv ~ 2 MeV. For example, with 
these approximations the decay amplitude for the GR model becomes A{GR) 
[8^/2/(27^)3]^eff («^^ - w''^^)ppEp. 

2 To perform the phase-space integrals, Pn = / (i£i(i£2£i£2(Q~^i "^2)", we neglect both 
the electron mass, me, and the Fermi functions, F{ei). The integrals then become: 
Pi = g^l, 260; P3 = <5Vl5, 120; P5 = Q"/83, 160; and P7 = Q^V308, 880. 

With these estimates, the effective equivalent GR coupling for which the (3f3^p decay 
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rate equals the f3f3^p^p decay rate is 



Qe^if^M ~ (0.008) g\sl4 




10-^ (10-^). 



(27) 



With eq.(22) (respectively (23)) in mind, we take gj^ ~ 0.2 (respectively Qj^ ~ 1), se ~ 0.1, 
Q ^ 2 MeV, m ~ 2 MeV (respectively m ~ 350 MeV) and pp ~ 60 MeV in obtaining 
this number. Although a coupling of 10~^ is probably hopeless to observe in present 
experiments (which can detect Qea ^ lO"'^), one at the level of 10~^ is close enough to 
observability to warrant a more detailed comparison of this type of model with the data. 

3.2) Decays Mediated by Virtual Scalars 

We next turn to models for which the PP^pip decay is mediated by the self couplings 
in the scalar sector, without the need for any extra heavy sterile leptons. We argue in 
this section that PP^pip decay in these models is always constrained by LEP data to be 
too small to be observed. To see how this works, consider an example of such a model in 
which the standard model is extended by the inclusion of: (i) a second SM Higgs doublet, 
X, which has non-zero lepton number L = — 1; (ii) a singlet Higgs boson field, 0, with 
lepton number L = +1, and {in) a weak isotriplet scalar field. A, with L = —2. These 
charge assignments permit the couplings considered in the previous section: the A Yukawa 
couplings contain a term of the form L L A, which couples it to the electron, and the Higgs 
potential contains the following terms which are relevant to PPip,p decay: 



Provided that the singlet field, (p, is sufficiently light, the diagram of Fig. (3) then 
induces PP^p^ decay. In this model the extra scalar fields do not have nonzero vev^s so that 
lepton number remains a good symmetry and the usual neutrinos remain massless, as in 
the SM. 

Neglecting E'^Wp in comparison to p^w^^, we find the approximate form for the 
amplitude of eq. (16): 



(28) 




(29) 
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where sq oc fii is the mixing angle which controls the strength of the mixing between the 
light mass eigenstate, (p, and the singlet field, 0. J^{pf, M^) is the function obtained by- 
performing the £ integration, which takes the values ^ ~ 1 for M,^ < pp and ~ {pp/M-^Y 
for ^ Pp. With these expressions we find: 

ge^siPM ^ (0.02) Jf-f hsl (^) :Fipp, M^). (30) 

As usual, the rate is largest if all of the exotic scalars have couplings that are as large 
as possible, and masses that are in the range of (10 — 100) MeV. This range is particularly 

dangerous for the phenomenology of the present model, however, since the new scalars 
couple to the photon and to the Z. Any such particles having masses < 50 GeV would 
contribute unacceptably to e+e~ annihilation, and so are ruled out. The effective coupling, 
(7efr, for scalar- mediated /3/3<^(^ decay rate is therefore suppressed at least by the factors 
(Q/Ma) {pp/M^Y ~ 10~^^, and so is much too small to be detectable. 

Since the suppression factor due to the large scalar masses is so devastatingly small, 
this rules out observable scalar- meditated |3(3^p^p decays in a much larger class of models 
than those we have directly considered. That is, one might entertain models for which 
scalar-mediated f3|3^p^p decay does not involve the W boson at all. After all. Fig. (3) works 
equally well if the longitudinal parts of the W lines are replaced by yet another exotic scalar 
which couples directly to quarks. In this case the nuclear matrix elements can be larger 
since they can involve quark operators other than the electroweak charged current. Also, 
the Yukawa coupling of such a scalar to quarks could easily be of order 10~^ — a value 
which is much larger than the corresponding value of 10~^ for the SM Higgs. Nevertheless, 
none of these enhancements can overcome the suppression due to large scalar masses. 

4. Model-Building: The Case of Spectral Index n — 7 

We now turn to a model which produces /3/3 decay with spectral index n = 7, i.e. 
a model of type HE. To do so we require a model in which lepton number is unbroken, 
and which contains a Goldstone boson, (/?, that carries lepton number +1. This can only 
happen if the theory has a nonabelian flavour symmetry which acts on leptons, and which 
contains ordinary lepton number as a generator. 

The simplest case is to take the lepton-number symmetry group to be G = SU (2) x 
Ul'{1), with G broken down to the Ul{1) of lepton number. This can be done by working 
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with a variation of the model of ref. [5]. We therefore add the following electroweak-singlet, 
left-handed fermions to the standard model: N = (^^), so and s_, which transform under 
the flavour symmetry, G, as: N ~ sq ~ (1,0) and s_ ~ (1,-1). If the unbroken 

lepton number is chosen to be Lg = T3 + L', where L' is the generator of [/^^/(l), then 
the subscripts of N+, Nq, sq and s_ give the corresponding particle's lepton charge. In 
order to implement the symmetry breaking pattern G 1/^(1), we also introduce the 
electroweak-singlet scalar field, $ = (^+) , which transforms under G as: $ ~ (2, . 

For this model the most general renormalizable G-invariant mass and Yukawa inter- 
actions involving the new fields are 

>^^m = - ylsoTiSo) + h.c., ^^^^ 
Cy = -\{s--i^L) H-g. (iV,7^6'_) e'^ + go (N^-,'r.so) e'^ + h.c. 

If the scalar fields acquire the vev^s {H^'^ = v and {(po) = u (which we take, for 
simplicity, to be real), then G breaks to ^7^(1) as required. If we also neglect mixings 
with and the following neutrino spectrum is produced. First, there is one massless 
Le = +1 state, z/g = UeCe — N+sg, where the mixing angle, 6*, satisfies: tan 6* = Xv/g-U. 
Then, there is a massive Lg = +1 Dirac state, Us — 7i-^+ +7jiS-, with = N+ce + i'eS0, 
and with mass Mg = y^A^w^ + g'^u'^. There are also two majorana L^, — states, i'±, 

which respectively have masses M± = I [yM2 + Ag^u'^ ± m] . These mass eigenstates are 
related to Nq and sq by: 



, . _ , . , . (32) 

f-j- / y Sip J \ Sq / 

with tan(2(^) = 2gQu/M. The factors of i in this mixing matrix are due to the chiral 
rotation needed to ensure that both mass eigenvalues, M±, are positive. 

The next step is to rotate fields to make the derivative couplings of the Goldstone 
bosons explicit (this is done in detail for a related model in ref. [5]). The appearance of 
the derivative coupling can be seen as follows: use exponential parameterization of the 
scalar fields that have non-zero vev so that the Goldstone field appears in the exponential. 
Then redefine the fermion fields so that the new fermion field is given by F' = Fe'^/'' (in 
this case F' is the heavy neutral lepton). After this redefinition, the Yukawa interaction 
becomes completely independent of the Goldstone field, which reappears with a derivative 
coupling due to the fermion kinetic energy term. The final form of the interaction given in 
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eq.(8) then emerges on diagonalization of the fermion mass terms with couphng matrices 
Xia and Yia given by Yia = 0, and 

_.S9C^ ^ _ S0s^ _ cec^ _ ces^ 

u " u u u 

With these expressions we may evaluate the PP,^,^ matrix element of eq. (10). It 
becomes: 



r d*. 
J (2^ 



2 slclMf 

(34) 



X 



'4 (^^ - ^e) (^2 + ^1 - i^Y 



4M+ 4m_ 



(£2 + M2 - ie) (£2 + M2 - ie) 



We quantify the size of this rate in a similar way as already done for the model of 
section 3.1, namely by means of estimating the magnitude of (/eff. We therefore choose 
values for masses and mixing angles that are compatible with the phenomenological bounds 
of section 5, but which also yield a maximal rate. Specifically we take 

Ms > 350 MeV; M+, M_ ~ a few MeV; ^ ~ 7.5 x 10"^; g- - 0.2. (35) 

The amplitude can then be written in the following approximate form: 



74 



X 



(36) 



(27r)4 (£2 _ ie) (£2 + m2 - ie)(£2 + - ie)' 



To see this keep in mind that in the limit of small 9 we have u ~ Comparing this with 
the amplitude of the GR model gives an equivalent of 



9- 



^eff ~ (0.0006) slclg- (j^) - 10"'"- (37) 

where the values of eq.(35) have been applied and we take the mixing angle to satisfy 
tan {(p) ~ l/-\/2 to maximize the rate {i.e. equivalently M+ ~ 2M_). Specifically we take 
M+ ~ ~ Q ~ 2 MeV, g- ~ 0.2, ~ 0.1 and Mg ~ 6pj. ~ 360 MeV. 
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In the above we have chosen 0.2 as the maximal value for as is naively required 
by kaon decay measurements (see section 5). It is probable that values as large as ^_ ~ 1 
would actually prove to be consistent with these measurements in a more careful treatment, 
however. This is because the scalar's derivative coupling provides an additional suppression 
to its contribution to kaon decays. If so, then the upper bound on g^fi can be relaxed to 
around ^feflf < 10~^^. 

It is clear that this rate is far too small to have any chance of being detected, largely 
due to the suppression by powers of Q that originate with the extremely soft electron 
spectrum. We consider this model therefore only as an example that such an exotic decay 
spectrum, with spectral index n = 7, can in principle exist. Given that this is also the 
kind of spectrum that is least well constrained experimentally, since cuts usually exclude 
the lowest-energy electrons which are usually the most contaminated by backgrounds, it 
would be very interesting to construct a model which predicts a detectable rate. We have 
so far been unable to find such a model. 

5. Phenomenological Constraints 

The previous sections have seen us use particular kinds of masses and couplings in 
order to maximize the decay rate. The present section is meant to justify the masses 

and couplings we have chosen. We do so by listing the main phenomenological constraints 
which the models we consider must confront. Since the models with the largest signals 
involve sterile neutrinos coupled to light scalars, we focus principally on the constraints on 
these. 

5.1) Laboratory Limits 

In most of the models considered in this paper the SM has been supplemented by a 
number of electroweak-singlet fermion and scalar fields. Of these, it is only the mixing 
of the sterile neutrinos with (in this case) the electron neutrino, u^, that is bounded 
experimentally. Since the sterile neutrinos in the models we consider dominantly decay 
into lighter neutrinos and scalars, all of whom are invisible to present day detectors, the 
only bounds which need be considered are those which do not assume a decay into a visible 
mode. The bounds that do apply can be classified into two subclasses according to whether 
they are due to low- or high-energy experiments [20]. 
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At low energies we have those experiments which examine the decays of pions [21] 
and kaons [22] at rest, and search for nonstandard contributions to the decay rate and the 
outgoing electron spectrum. These experiments constrain the masses of sterile neutrinos 
in a model-independent way from a few MeV almost up to ~ 100 MeV for pion decays, 
and up to ~ 350 MeV for kaon decays. Mixing angles, Uei, with the electron neutrino, 
i^e, are strongly constrained for sterile neutrinos with masses in this range. For example, 
upper bounds on \Uei\ can be as low as 7 x 10"'^ for a mass of 50 MeV [21]. 

These same experiments, as well as measurements of the muon lifetime, also limit the 
possibility for the existence of weak decays into very light scalars in addition to neutrinos 
[23]. Specifically these limits come from precision measurements of the Michel parameter, 
p, in muon decay, and from using the measured decay spectra of the outgoing leptons to 
constrain decays such as K — > £N(p or tt — > £N(p. The failure to observe such decays 
can be expressed as an upper bound on a hypothetical dimensionless Yukawa coupling, 
^eff, of an effective — N — ip interaction. In terms of this coupling, the current bounds 
are respectively ^eff ^ 5.7 x 10~^ (Muon decay), < 1.6 x 10~^ (Pion decay) and < 1.3 x 
10~^ (Kaon decay) [23].^ Of course, these bounds assume that the kinematics permit the 
emission of N^p in these decays. For the models of interest here, the bound therefore applies 
if the Le = sterile fermions, Na, are sufficiently light compared to the tt, K and/or |U, as 
is certainly the case for mjv„ ~ 1 MeV. For mjv^ in this mass range, the effective coupling 
that is then bounded in this analysis turns out to be ^eff ~ sgceg^. As a result, keeping 
in mind the bound on 9 discussed below, g+ must be smaller than 0.18. 

It is clear however that provided Na is too heavy to be produced in these kinds of 
decays (i.e. mjv„ ^ 350 MeV) then the bound described above no longer applies, and g+ 
can be of order unity. We are led to consider m^^^ ~ 1 MeV because of the specific way we 
choose to avoid conflict with nucleosynthesis (see section 5.2). If another way to evade this 
bound should be found, or if we simply accept the comparatively large number of effective 
degrees of freedom contributed at this epoch, then g+ can indeed be as large as g+ ^ 1. 

For masses heavier than 350 MeV the best bounds on sterile neutrinos usually come 
from beam dump experiments. However this type of experiment relies crucially on the 
assumption that the isosinglet leptons decay via their charged current weak interactions. 



These bounds can be even more stringent if additional model-dependent assumptions for the couplings 
involved are made. 
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The same is also true for bounds coming from searches on the Z-pole for exotic decays of 
the Z boson, such as Z — > Nu — > W*ev. (Here AT denotes the isosinglet state, and W* 
denotes a virtual W boson.) However, since such decays do not dominate in the models we 
consider — the sterile states instead prefer to decay into light scalars — these experiments 
need not be considered further here. 

The best limits, for the mass range above ~ 350 MeV, therefore come from the re- 
duction of the effective charged- current coupling of the electron, which is suppressed by 
the cosine of the mixing angle between i^e cind the heavy sterile state. This reduction 
potentially affects electroweak precision experiments by influencing the measured value of 
Fermi's constant, G^., that is inferred from muon decay. It also shows up as a failure of 
lepton universality in low-energy weak decays. This leads to the bound \U^i\ < 7.5 x 10~^ 
(2cr) for masses above ~ 350 MeV [24]. Another bound comes from the reduction of the 
neutral current coupling of the electron-neutrino, leading to a decreased number of light 
neutrino species as deduced from Z's invisible width. However this bound is somewhat 
weaker than the aforementioned universality bound, and applies only to sterile neutrinos 
that are heavier than about 90 GeV, hence it is of no particular importance in our case. 

These bounds have immediate implications for the models considered here, such as 
those discussed in sections 3.1 and 4. Since these models have only one sterile state, 
the heavy Dirac neutrino Us, that mixes with u^, we need only constrain the masses and 
mixings of this one particle. To maximize the /3/3<^<^ rate we prefer a mixing angle, 9, that 
is as large as possible, and so we choose Vg to be heavy enough to evade the bounds from 
pion and kaon decay — Mg > 350 MeV — and we take 9 to saturate the universality 
bound — ^ < 7.5 X 10-^. 

5.2) Nucleosynthesis 

Standard big-bang cosmology very successfully explains the abundances of light ele- 
ments produced during the nucleosynthesis epoch [25]. In fact, this success leaves little 
room for the existence of new species of particles at this critical time. Specifically if there 
were new particles present at a temperature of ~ 0.1 — 2 MeV their additional degrees of 
freedom would increase the expansion rate of the universe leading to an earlier freeze-out 
of neutron-proton converting interactions and thus to a larger amount of neutrons that 
eventually would be cooked into "^He. The bound on the number of additional degrees 
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of freedom is conventionally stated in terms of the maximal number of additional light 
neutrino species allowed at nucleosynthesis 5Nj^ [26]: 

5N^ < 0.4 (38) 

There are two simple ways to satisfy this bound when additional particle types are 
present in a theory. For particles that are heavier than about 1 MeV, the most straight- 
forward way is to have them not be present at all during nucleosynthesis. This can be 
achieved by either having them decouple and subsequently decay sufficiently early — keep- 
ing in mind the possible contributions of the decay products — or by having them remain 
in thermal equilibrium as they become nonrelativistic, so that their abundance becomes 
suppressed by the Boltzmann factor, e""*/^. For particles which are much lighter than 1 
MeV, on the other hand, the contribution to 5N^ is strongly suppressed if they decouple 
before the QCD phase transition at T ~ 200 MeV. In this case they become sufiiciently 
diluted compared to ordinary particles afterwards due to the reheating of the photon bath 
during this phase transition. 

Unfortunately, these methods alone are not enough to avoid conflict of the models 
considered in this paper with nucleosynthesis. This is because the very light scalars, 
that are present in all of our models, couple through dimcnsionlcss couplings to ordinary 
neutrinos, such as Ug. As a result, they remain in equilibrium with these neutrinos down 
to temperatures that are well below 1 MeV. This makes a conflict with nucleosynthesis 
generic for models which can produce PP^, and PP^p^ decays. 

Rather than therefore considering these models to be inevitably ruled out, we demon- 
strate here an existence proof that the nucleosynthesis bound can be avoided, given suitable 
masses for the sterile neutrinos. We propose this as an example of the kinds of informa- 
tion that could be inferred from cosmological bounds, should one of the exotic decays 
we describe ever be detected experimentally. The loophole we describe follows its similar 
application in ref. [18], and is based on an observation of ref. [27]. 

Suppose, then, that a heavy neutrino state either decays [27] or annihilates [18], into 
particles that are in equilibrium with ordinary neutrinos, right after these neutrinos freeze 
out of chemical equilibrium with the photon bath {i.e. T^, ~ 2.3 MeV.) Provided this hap- 
pens before the neutron-proton ratio freezes out (at T^/p ~ 0.7 MeV), then the abundance 
of neutrinos will be increased relative to the standard cosmological scenario. This over- 
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abundance of neutrinos will act to suppress the neutron abundance at freezeout, and so 
decreases the predicted production of ^iJe. As a result it contributes to nucleosynthesis 
as a negative contribution to SN^, [27] . 

Consider, therefore, a single majorana neutrino with a mass in the critical region [i.e. 
~ 1 — 2 MeV) which is still in equilibrium with Us species of real scalars in addition to 
the Up — 2) species of ordinary neutrinos, when it becomes nonrelativistic. Such a particle 
will annihilate out right during the critical epoch, thereby heating the neutrino-scalar bath 
above that of the photon bath. Then the magnitude of SN^ that results due to the above 
mechanism turns out to be [27], [18] SN^, = — 4.6 5n, where 

Here n/uo and E/Eo denote the ratios of the electron- neutrino number densities and 
energies after and before annihilation. The second expression is obtained from the first 
one using entropy conservation. 

We can now see how the models described in previous sections can accommodate the 
nucleosynthesis bound. Take the model of section 3.1 as a representative case. This model 
contains a massive Dirac state, z/^, which mixes appreciably with the electron neutrino, 
and so which must be heavier than ~ 350 MeV to avoid the laboratory bounds described 
earlier. It also contains a massive majorana neutrino, A^O: as well as a very light complex 
scalar. 

Since the sterile neutrino is kept in thermal equilibrium with the photon bath by 
neutrino scattering via scalar (0) exchange, it remains in equilibrium long after the tem- 
perature has dropped below its mass. Its abundance is therefore sufficiently suppressed to 
ensure that it never plays a role at nucleosynthesis. The same would be true for the ma- 
jorana neutrino if it were heavy, however this would leave the light scalars in equilibrium 
at nucleosynthesis, contributing SN,^ = |. The contribution of these scalars can be very 
efficiently cancelled, however, if Nq should have a mass of around 1 MeV, since in this 
case it annihilates right in time to reheat the neutrino sector.^ In this case we can employ 
eq. (39) with = 3 and ns = 2, yielding 5n = 0.43 and so a total shift in the effective 
number of neutrino species of (5iV^)tot = f- 4.6(0.43) = -0.86. Since (5iV^)tot rises to f as 

No remains in equilibrium with Ue through scalar exchange. 
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A^o's mass increases above Tj, ~ 2.3 MeV, it is possible to find an efi^ective number of fight 
neutrinos which is compatible with observations. To summarize this paragraph: Confiict 
with nucleosynthesis can be avoided if A^o has a mass in the range of a few MeV. 

Following similar lines of reasoning one can also make the model of section 4 satisfy the 
nucleosynthesis bound. In this case the annihilation of just one majorana neutrino during 
the critical epoch is not enough though because we have twice as many light scalar degrees 
of freedom as in the model discussed above. Assuming however that the two majorana 
mass eigenstates u± are both not much heavier than a few MeV solves the problem. To 
make this point explicit we just state that if both and annihilated out right in the 
critical epoch then this would result in (5A^z/)tot = — 1- 

5.3) Other Constraints from Astrophysics and Cosmology 

Other astrophysical bounds can constrain the properties of very weakly interacting 
particles, such as the sterile neutrinos and scalars we are entertaining. These turn out to 
not significantly constrain the models we consider. We illustrate this point in this section 
by examining these bounds for the model considered in section 3.1. 

• (1) Stellar Evolution: Sterile particles produced in the core of stars can disrupt our 
understanding of stellar evolution since they can provide an extremely efficient additional 
cooling mechanism for a star. For the model of section 3.1, the sterile neutrinos Vg and A^o 
are simply too heavy to be produced in significant numbers inside the core of a star, since 
their masses are respectively taken to be ~ 350 MeV and ~ 1 MeV. They are therefore 
not constrained by stellar evolution. The same conclusion also holds for the light scalars 
of this model, although for different reasons. In this case the scalars are light enough 
to be produced, but they only couple appreciably to neutrinos. But since temperatures 
and densities inside a star are not sufficiently high to produce a thermal or degenerate 
population of neutrinos, scalars are not produced in significant quantities from the interior 
of the star. Direct couplings with other particles, such as electrons, do arise due to loop 
effects, but these are too small to be significant.^ 

• (2) Supernovae: Populations of neutrinos are maintained, however, within the cores 
of supernovae, due to the much higher temperatures (~ 50 — 100 MeV) and densities 

^ See, for example, the second Ref. of [5], where this issue is treated in more detail in a similar context. 
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that are found there. Both the scalar, and the sterile majorana neutrino, A^O) cire 
therefore sufficiently light to be produced in significant numbers. However since both of 
these particles are also in thermal equilibrium with the ordinary SM neutrinos, they are 
trapped in the core of the supernova and so do not lead to premature cooling. 

• (3) 'Present-Day' Cosmology: One last issue is the effect of (p on structure formation, 
and on the energy density of the universe observed today. Since (p would remain in thermal 
equilibrium with the electron neutrinos (through exchange of virtual A^o) until a tempera- 
ture in the region of 1 eV if it was light, it inevitably annihilates out to the neutrino sector 
once it becomes nonrelativistic. It will therefore never dominate the energy density of the 
universe and hence also not affect structure formation nor the age of the universe in any 
significant way. This in turn means that (p can have any mass below 1 MeV. 

6. Conclusions 

In this paper we investigate a novel class of neutrinoless double-beta decays, in which 
two (rather than one or no) scalars are emitted in addition to the two observed electrons. 

We do so to provide sample models of the two remaining classes of spectra for PP decays, 
which had not been hitherto explored. We regard it as useful to construct representative 
models for each class, since such models are required in order to determine the implications 
of other experiments, at different energies, for (3(3 decay. 

We have obtained the following results: 

• 1: We have found that two-scalar decays always lead to a (3(3 electron sum-energy 
spectrum whose spectral index is either n = 3 or n — 7. The n = 3 spectrum is intermediate 
between the SM shape (for which n = 5), and the traditional scalar-emitting shape (n = 1), 
such as is found in the Gelmini-Roncadelli model. The n = 3 spectrum is identical to a 
recently-discovered class of models [5] for which only a single scalar is emitted during the 
decay. 

• 2: The n = 7 spectrum is completely new, and is softer than the observed SM (3(32u 
spectrum. As a result, it is more difficult to experimentally extract, since most backgrounds 
tend to dominantly produce low-energy electrons. It is intriguing that such an exotic 
spectrum could in principle be hiding undetected even now in the present /3/?2i/ data! 
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• 3: Viable models can be constructed which predict an n = 3 (3(3 decay rate which is 
reasonably close to the present experimental limit. In the viable models with detectable 
rates the decay is mediated by the exchange of a number of sterile neutrinos, whose masses 
must lie in the (1-100) MeV range. Their couplings must also be as large as is permitted. 
The requirement of such large couplings to obtain observable signals suggests that an 
observed I3[3^p^p signal should be accompanied by the appearance of new physics in other 
experiments, such as in violations of lepton universality. 

A potential problem with this model is that it predicts an excess number of degrees 
of freedom at nucleosynthesis that is equivalent to | of a neutrino, in comparison with the 
presently-quoted bound of 0.4. A way of evading this nucleosynthesis bound is described 
in section 5, but the required couplings and masses for this scenario would predict a 
rate which is two orders of magnitude smaller, and so out of reach of current experiments. 
We do not regard this as sufficiently worrying to preclude searching for this decay in 
current experiments, since an experimental signal would likely prompt more imaginative 
approaches to understanding nucleosynthesis in these models. 

• 4' Because of its soft spectrum, we were unable to construct viable models of the n = 7 
P(3^^ decay, for which the predicted rate is large enough to be detectable. Thus, although 
the existence of an n = 7 decay mode is experimentally tantalizing, it is unlikely to be 
discovered in decay experiments for the foreseeable future. 

• 5: For the purposes of comparison we also examined a class of models for which PP^^ 
arises due to the mutual interactions amongst new particles in the scalar sector, rather 
than from new sterile leptons. Although the decay rate that is predicted for these models 
is also too small to be seen, the decays nevertheless exhibit novel features. One such is a 
nonstandard angular distribution of the two electrons as a function of the opening angle 
between them, similar to what arises in models with right-handed currents. 

All of these results serve to underline the theoretical understanding which has emerged 
over the past few years. Scalar-emitting [3/3 decays can reasonably be expected to be 
observed, although the properties of the phenomenologically viable models which can do 
so are very different than what would be expected based on intuition that is based on the 
original GR-type models which originally motivated these experiments. 
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8. Figure Captions 



(1) The (3(3 spectrum as a function of the two electrons' total kinetic energy for various 
choices of the 'spectral index' n. n = 1 corresponds to the dotted line, n = 3 is the 
dashed line, n = 5 is the solid line and n = 7 is the dash-dotted line. All four curves 
have been arbitrarily assigned the same maximal value for purposes of comparison. 

(2) The Feynman graph representing the leading order contribution to the two-scalar 
/3/3<^<^ decay through sterile lepton exchange. 

(3) The Feynman graph contributing to (3(3ip^p decay purely due to the self-interactions of 
the Higgs sector. 
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